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Abstract. The millimeter and sub-millimeter waves have 
been attracting a lot of attention recently in the cloud re- 
mote sensing community. This is largely because of their 
potential use in measuring cirrus cloud parameters with air- 
borne or space-borne radiometers. In this study, we examine 
the possibility of using polarization measurements in this fre- 
quency range to get information on the microphysical prop- 
erties of cirrus clouds. By using a simple radiative transfer 
model, we calculated the brightness temperature differences 
at the vertical and horizontal polarization channels for the 
following seven frequencies: 90, 157, 220, 340, 463, 683, 
and 874 GHz. The ice crystals in cirrus clouds are modeled 
with nearly spherical particles, circular cylinder, and circu- 
lar plate, as well as with mixtures of these types. We found 
that the polarization difference signal shows a unique “reso- 
nance” feature with the change of ice particle characteristic 
size: it has a strong response only in a certain range of ice 
particle size, beyond that range it approaches zero. The size 
range where this resonance happens depends to a large extent 
on particle shape and aspect ratio, but to a much less extent 
on particle orientation. This resonance feature appears even 
when ice clouds are composed of a mixture of ice crystals in 
different shapes, although the magnitude and the position of 
the resonance peak may change, depending on how the mix- 
ture is made. Oriented particles generally show larger po- 
larization difference than randomly oriented ones, and plates 
have larger polarization difference than cylinders. However, 
the state of particle orientation has a significantly stronger 
effect on the polarization difference than the particle shape 
(cylinder or plate). This makes it difficult to distinguish par- 
ticle shapes using millimeter and sub-millimeter radiometric 
measurements, if there is no information available on parti- 
cle orientations. However, if the state of particle shape mix- 
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ture can be predetermined by other approaches, polarization 
measurements can help to determine ice particle characteris- 
tic size and orientation. This information, in turn, will benefit 
our retrieval of the ice water path of cirrus clouds. 


1 Introduction 


In the last two decades many studies have been carried out for 
the remote sensing of cirrus clouds from ground, airplane, 
and satellites (Lynch et al., 2002). Both active and passive 
approaches have been tested. Most of them are exploiting 
the advantages of the visible, infrared, and microwave range 
of the electromagnetic spectrum. Recently, the millime- 
ter (mm) and sub-millimeter (sub-mm) wave range, which 
covers frequencies from about 150 to about 1000 GHz with 
corresponding wavelength ranging from 2.0 to 0.3 mm, has 
been explored. Both theoretical simulations and experimen- 
tal studies have been conducted. The results show that this 
spectral range could be very promising in measuring cloud 
parameters such as the ice water path (IWP) and the ice parti- 
cle characteristic size (Evans et al., 1998; Wang et al., 2001). 

The mm and sub-mm range shows some important ad- 
vantages compared with the other commonly used spectral 
ranges. First of all, the wavelengths in this range are com- 
parable with or larger than the sizes of most ice particles. 
Secondly, the imaginary part of the refractive index of pris- 
tine ice in the mm and sub-mm spectrum is much smaller 
than the real part of the refractive index. These two factors 
combined cause the interaction between ice particles and the 
atmospheric radiation to be dominated by the scattering ef- 
fect, which is almost independent of the ambient temperature 
of the cloud and mainly related to the volume (or mass) of the 
ice particles. Furthermore, because of the strong absorption 
of atmospheric water vapor, ground surface is almost unde- 
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tectable for a space-borne mm and sub-mm radiometer. Fi- 
nally, since cirrus clouds are mostly thin, the brightness tem- 
perature depressions induced by cloud ice particle scattering 
are, in most cases, proportional to the cloud IWP, at the mm 
and sub-mm frequencies (Liu and Curry, 2000). This will 
significantly simplify the retrieval of cloud parameters from 
radiometric measurements. 

However, until now most studies are concerned with the 
information content of radiometric intensity signals (Evans 
et al., 1998; Wang et al., 2001). Frankly speaking, there is a 
general lack of studies on the potential of polarization mea- 
surements in the mm and sub-mm wave band. In this paper, 
we will treat this topic in some detail. As mentioned ear- 
lier, the ground surface is mostly undetectable for a space- 
borne mm and sub-mm radiometer, except for the two lowest 
frequencies at 90 and 157 GHz considered in this paper, at 
which the ground surface contribution can be significant in 
dry atmospheric situations. In our later simulations, we ig- 
nore the surface contribution at these two frequencies and 
at other frequencies as well. Under this condition, the po- 
larization signal is exclusively induced by the scattering of 
frozen hydrometeors, i.e. ice crystals to the unpolarized at- 
mospheric radiations. By using a simple physical model, we 
demonstrate in Sect. 2 the advantages of polarization mea- 
surements. In Sect. 3 we apply this model to cirrus clouds 
in order to examine factors which influence the polarization 
signature measured by a space-borne mm and sub-mm ra- 
diometer. These factors include particle shape, particle ori- 
entation, and particle size. And finally in Sect. 4, we present 
some discussions and conclude the paper. 


2 The physical model 


The polarization signature is described by the brightness tem- 
perature difference measured by a radiometer at two orthogo- 
nally polarized channels, i.e. the vertical and horizontal chan- 
nels. Highly accurate polarization difference measurements 
at the microwave band can be achieved with today’s tech- 
nology (Hollinger et al., 1990). The usage of this signature 
has been shown to be very effective, e.g. in determining sur- 
face features as sea ice concentration (Kern and Heygster, 
2001) or atmospheric parameters such as rain rate over land 
(Spencer et al., 1989) and cloud liquid water over sea ice 
(Miao et al., 2000). To analyze the information content of 
the polarization signature of cirrus clouds, we use the simple 
physical model shown in Fig. 1. 

In this simple model, the atmosphere is assumed to be in 
a plane-parallel structure containing only one single layer of 
ice cloud. The atmosphere has three layers, which are as- 
sumed to be homogeneous and isothermal. Only the middle 
layer contains ice particles. The brightness temperature ob- 
served by a space-borne radiometer at the polarization state 6 
is denoted by Tpg and £ can be v or h, respectively for verti- 
cal or horizontal polarization. Tpg can be expressed explicitly 
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with the physical parameters given for each layer. The reflec- 
tion and refraction effects at layer boundaries are neglected. 
Within the cloud layer only single scattering effect is taken 
into account for the interaction between ice particles and at- 
mospheric radiation. Under these conditions, Tpg consists of 
the following four parts: the first part, denoted by Tyg 1, is 
the emission from gases in the uppermost layer. The second 
part, Tpg,2, is the emission from the layer containing a cloud, 
which is attenuated by gaseous absorption in the uppermost 
layer. The third part, Tpg,3, represents the radiation coming 
from the lowermost layer, which is attenuated by the overly- 
ing two layers through both absorption and scattering. And 
finally the fourth part, Tpg,4, describes the contribution due 
to ice particle scattering in the cloud layer. This part is also 
attenuated by the upper layer absorptions. Summing up all 
these parts gives 


Tog = Top,1 + Tog,2 + Top,3 + Top,4 (1) 


The terms on the right-hand side of Eq. (1) can be explicitly 
expressed as 


Tog 1 = M1 —e ™! sec 0) és, 
Tp.2 = Ta (1 — e122 889) etal Sec G 
Tog. 3 = DA — e sec 0) o~ (Tal + tart ts2,p) sec 8 (4) 
and 

Tpp,4 = Tyge “4! sec 0 6 


Here Tai (i € 1, 2, 3) is the partial opacity of layer i due to at- 
mospheric absorptions. Theoretically, t,2 includes the effect 
of cloud ice particle absorptions and consequently it depends 
also on the polarization state. However, since ice particles 
absorb little mm and sub-mm wave energy, we hereby ne- 
glect their contribution in Ta2. The attenuation caused by ice 
particle scattering is represented by Ts2,g, which depends on 
the polarization state 6. The physical temperature of each 
layer is denoted by T; with i € 1,2,3. Additionally, Tsg 
in Eq. (5) represents the scattering contribution of ice parti- 
cle scattering, which was made by radiations coming from 
other directions and scattered toward the radiometer observ- 
ing direction. This term is directly related to ice particle mi- 
crophysical properties as particle number, size, shape, and 
orientation. 

Since the first two parts of Tpg, i.e. terms Tpg,ı and Tyg 2, 
do not depend on the polarization state, these two terms can 
be canceled by taking differences of the brightness tempera- 
tures measured at two orthogonally polarized channels, 


Q = Thy — Ton. (6) 


In this way, the polarization difference Q will not depend on 
atmospheric emissions from the cloud and from the upper- 
most layer, although the uppermost layer still have influences 
on Q through the absorption attenuation. This is because of 


www.atmos-chem-phys.org/acp/3/39/ 


J. Miao et al.: Polarization measurements of cirrus clouds 


41 


= 
—_—_. Ə $ 
E 
P y JN: 
Tit 4 Thar oy OTs Tope 


Fig. 1. The geometry of the problem. 
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By defining Tpz as 


, which is contained in terms Tyg 3 


Tor = T3(1 — e 73 °9), (7) 
can be explicitly written as 
p y 
Q Z The Tait tar) sec 0 (e7Ts2v secO eT Ts2h sect 


+(Tsy — Tsne ™ basse (8) 


Note that Tpz approaches to the layer air temperature T3 
as the layer opacity t,3 approaches to infinity. The right- 
hand side of the above equation shows that Q consists of two 
terms, which describe the contributions from two different 
sources. The first term represents the direct contribution of 
the radiation coming from the lowermost layer and traveling 
through the cloud and the uppermost layer up to the space- 
borne radiometer. Since the radiation from the lowermost 
layer is unpolarized (note, the ground surface polarization 
signal is ignored), the polarization difference is exclusively 
induced by the different scattering efficiencies of ice particles 
at the two different polarization states. The second term, in 
contrast, expresses the contributions of ice particle scattering 
of the radiation from all directions. The first term will be zero 
if ice particles are spherical or randomly oriented, because, 
in this case, ice particle scattering efficiencies are same at the 
two orthogonal polarization states, i.e. T;2,y is equal to Ts2,}. 
In general, however, the second term is nonzero. 

Suppose that we have two different frequencies, fı and 
jo, at which the opacities due to atmospheric absorption are 
same, i.e. Tai( f1) is equal to Tai( f2), for all three layers 
(i € 1,2,3), then the ratio of the two polarization differ- 
ences gu will be free of the screening effect (due to both 
emission and attenuation) of the upper layer. By referring 
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to Eq. (8), we can understand that this is because the two 
terms on the righthand side of Eq. (8) contain a common fac- 
tor e— 18°C? Furthermore, since Q is proportional to TpL, 
which represents the intensity of radiation coming from the 
lowermost layer, the ratio a : is, to some extent, also in- 
dependent of the atmospheric conditions below the cloud. 
This will simplify our consequent cloud parameter retrievals. 
Based on this idea, a recently proposed concept of a space- 
borne mm and sub-mm radiometer, which is dedicated to 
cirrus cloud measurements, has included some well-matched 


frequencies (for details see Miao et al. (2002)). 


When the cloud is optically thin, we expect that the po- 
larization difference Q will increase linearly with the in- 
crease of cloud optical thickness (ice particle size and/or 
cloud IWP). If the cloud optical thickness increases continu- 
ally, Q will reaches its saturation due to multiple scattering 
at a certain amount of cloud optical thickness. How large 
this amount is depends on the observing frequency. There- 
fore, simultaneous measurements at two or more frequencies 
will enable us to cover a large range of cloud particle sizes 
and optical depths. On the other hand, for thin cirrus clouds 
where the polarization difference Q has a linear relationship 
with the cloud IWP, the evaluation of the ratio of polarization 
differences, Ofi) , can help us, at least to some extent, to get 
rid of the influence of IWP from the measurements and in 
such a way to further isolate the information of ice particle 
microphysical properties. This microphysical property infor- 
mation, when it is available, will help us in the retrieval of 
cloud IWP. 
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3 The information content of polarization measure- 
ments 


3.1 The radiative transfer model 


From the physical analysis in the preceding section, we know 
that the uppermost layer plays a minor role in the polariza- 
tion signature. Therefore we neglect the effect of the upper- 
most layer in the following discussions. By further assuming 
the cloud layer to be thin and omitting the absorption in this 
layer, we get a simpler form of Eq. (8) 


Q = Qs + Th (e729 808 — eB 20h SEO), (9) 


where Qs = Tsy — Tsn. From Mishchenko et al. (2000), we 
get Qs in the single scattering approximation, 


Q;(n) = No f (Zoi (n,n')) Tin nd. (10) 


Here, n and n’ denote the directions of the satellite obser- 
vation and the incident radiation, respectively. (Z2,(n, n’)) 
is the element (2,1) of the ice particle scattering matrix 
(Z(n, n’)), which is ensemble averaged according to ice par- 
ticle size, shape, and orientation distributions. No is the num- 
ber density of ice particles in the cloud times the cloud geo- 
metrical thickness and it is in a dimension of 1/m?. Jj, (n^) 
represents the intensity (in Kelvins) of radiation incident on 
the particle from a direction n’. 

To make the simulation more close to reality, we used the 
temperature and humidity profiles (with a surface relative hu- 
midity of 80%) of the standard US atmosphere (1976) (Ulaby 
et al., 1981), instead of the isothermal and homogeneous at- 
mosphere assumed in the physical model presented in the 
preceding section. The ground surface temperature is set to 
be equal to the surface air temperature and the surface emis- 
sivity is assumed to be 1. There is no polarization difference 
signal from the ground surface assumed. The exact value of 
the surface emissivity is relatively arbitrary, since the ground 
surface is undetectable for most frequencies considered in 
this study. The thin cirrus cloud layer is placed at the alti- 
tude of 10km. The satellite observation angle 6 is fixed at 
54°, which is mostly selected by low orbit satellites. Tpz and 
Tin(n’) are calculated using the radiative transfer model de- 
scribed by Gasiewski (1993). Since the atmosphere is taken 
to be plane-parallel, the intensity of the incident radiation, 
Tin(n’'), depends only on the incidence angle, which is the 
angle formed by the vertical direction of the altitude and the 
incident vector n’. At the altitude of 10km, the intensity of 
the radiation coming from the lower atmosphere is >230 K, 
which is much stronger than that from the upper atmosphere 
(a few tens of Kelvin) at all the frequencies considered in this 


paper. 
3.2 The microphysical model 


In situ measurements have shown that cirrus clouds are com- 
posed of various non-spherical ice particles (McFarquhar and 
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Heymsfield, 1996). For our simulations we selected three 
types of ice particles to represent them: the nearly spheri- 
cal particles, the circular cylinders, and the circular plates, 
for the ease of calculating the single scattering properties. 
All ice particles are assumed to be composed of pristine ice, 
which has a mass density of 0.92 gcm7?. In describing ice 
particle size distributions, we used the maximum dimension 
to represent the particle size, which is the length for cylin- 
ders and the diameter for plates. For nearly spherical parti- 
cles, the diameter of the mass equivalent sphere is used to 
define the particle size. In this paper, the maximum particle 
size is taken to be 2.0mm and the minimum is 20 wm. The 
size distribution is evaluated discretely for particle size from 
20 um to 2.0mm with a fixed step of 20 um. This range 
of particle size is reasonable, because ice particles smaller 
than 20 um are not significant enough for the mm and sub- 
mm wave scattering due to the very small volume fraction in 
clouds, although they may exist in a large number in realistic 
cirrus clouds (Lin et al., 1998). We assume the particle size 
distribution to be 
a@+3.67 D 


N(D) =aD%e~ Dn ? , (11) 


where a is a normalization constant and «œ is the width pa- 
rameter usually ranging from 0 to 2 (Evans et al., 1998). 
Parameter Dm is the median of the distribution of D to the 
third power, which separates the ice particles into two parts 
with equal mass, i.e. the total mass of ice particles with sizes 
smaller than Dm is equal to the total mass of ice particles 
with sizes larger than Dm. Since this distribution is discretely 
evaluated and also truncated at 2.0mm of D, there will be a 
mass loss for the part with sizes larger than Dm. This mass 
loss will increase with Dm. Our calculations show that, for 
a = 1, the mass loss is 1% of the total mass when D,, equals 
0.8mm. This loss increases to ~4% when Dm increases to 
1.0mm. The width parameter œ was found to play a mi- 
nor role in particle scattering properties, hence our following 
discussions will be limited to the case of a = 1 (Lemke and 
Quante, 1999). 

To calculate the single scattering parameters of ice parti- 
cles, i.e. parameters Ts2,8 (= No(Csca,p)) and (Z21(n, n’)), 
the T-Matrix method and the Discrete Dipole Approxima- 
tion (DDA) method are used (Mishchenko and Travis, 1998; 
Draine and Flatau, 1998). (Csca,g) is the ensemble-averaged 
scattering cross section of a single ice particle. Ensemble- 
averaging is done according to the given size, shape and ori- 
entation distributions. Because the T-Matrix method has a 
high computational efficiency, it is used in this study for most 
calculations of the single scattering properties. However, 
there are some cases when, perhaps because of the too large 
particle size and/or the aspect ratio, the T-Matrix code fails 
to converge. Then, the DDA method will be used to com- 
plement the T-Matrix method because it is more robust. The 
single scattering properties of the nearly spherical particles 
are calculated by averaging those from spheroid (both pro- 
late and oblate), circular cylinders and circular plates, which 
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Table 1. The refractive index of pristine ice at a temperature of 
250 K for the used mm and sub-mm frequencies 


f(GHz) A(mm) Real Part Imaginary Part 
90 3.333 1.7823 0.0026 
157 1.910 1.7818 0.0039 
220 1.364 1.7815 0.0049 
340 0.882 1.7815 0.0067 
463 0.648 1.7820 0.0084 
683 0.439 1.7846 0.0111 
874 0.343 1.7895 0.0134 


have aspect ratios (7) of less than 1.5. The aspect ratio 7 
denotes the ratio of the maximum to the minimum dimen- 
sion of a non-spherical particle. In situ observations have 
shown that ice particle aspect ratios can change from 1 to 
7 for cylinders and from 5 to 30 for plates (Auer and Veal, 
1970). Following this information, our simulations are car- 
ried out for circular cylinders with 7 from 1 to 7 with a step 
of 1 and for circular plates with discrete n of 3, 5, 10, 20, 
and 30. The frequencies used are at 90, 157, 220, 340, 463, 
683, and 874 GHz. Among these frequencies some (90, 157, 
220, and 340 GHz) have been already tested in space-borne 
or air-borne sensors and some (463, 683, and 874 GHz) are 
recently found to be well matched in the atmospheric absorp- 
tion comparisons and therefore they have been adopted by a 
new satellite mm/sub-mm radiometer concept (Miao et al., 
2002). The refractive index of pristine ice at the used fre- 
quencies is given in Table 1, which is taken from the data set 
given by Warren (1984) and for a specific ambient tempera- 
ture of 250 K. 


3.3 The information content of the polarization signature 
3.3.1 Cirrus clouds with randomly oriented ice particles 


In this section, we consider ice clouds composed of randomly 
oriented ice particles. Although aerodynamic and gravita- 
tional forces tend to cause ice needles and plates, especially 
those with large size and aspect ratio, to have a preferred hor- 
izontal orientation for their longer axes (Brussaard and Wat- 
son, 1995), the orientation angle distribution for the whole 
ensemble of ice crystals varies from case to case (Chepfer et 
al., 1999). Therefore, random orientation is often a good ap- 
proximation for ice particle orientations in cirrus clouds and 
it has been widely used in the retrievals of cirrus cloud pa- 
rameters such as particle shape and IWP (Baran et al., 2001; 
Chepfer et al., 2001). For randomly oriented ice particles, 
Ts2,y is equal to Ts2,,, Which results in a vanishing second 
term on the righthand side of Eq. (9). Therefore, only the 
first term Qs is needed to evaluate for the polarization sig- 
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874 Nearly Spherical Particles 


Polarization Difference Q (K) 


(0) 200 400 600 800 1000 
Diameter of Median Mass Equivalent Sphere Dm (um) 


Fig. 2. Simulated polarization differences for cirrus clouds com- 
posed of randomly oriented, nearly spherical ice particles. Results 
are shown for seven frequencies, which are in GHz and denoted 
above each curve. Detailed explanations are in the text. 


nature. Qs is zero for nadir-looking geometry and generally 
non-zero for non-nadir looking geometry. This is directly re- 
lated to the symmetry or asymmetry distribution of the inci- 
dent radiation intensity relative to the observation direction. 


In Fig. 2 simulation results are shown for the case of nearly 
spherical particles. Please keep in mind that the satellite ob- 
servation angle is 54°. The polarization difference Q is de- 
picted against the ice particle characteristic size, which is 
represented by the diameter of the median mass equivalent 
sphere. Seven curves correspond to results at the seven fre- 
quencies which were mentioned in the preceeding section. 
The cirrus cloud is assumed to have an integrated ice water 
path of 60 gm~*, which becomes 100 gm~? along the slant 
satellite observation line with a zenith angle 0 of 54°. Since 
only single scattering effect is considered in doing radiative 
transfer within the cloud, the simulation results given here 
can only be quantitatively used for thin clouds. Therefore, 
the values of the polarization difference Q given in Fig. 2 and 
in the following Figs. 3, 4, and 6 should be taken as the value 
per 60 gm~? of cloud IWP, when the cloud IWP is other than 
60 gm7?. 

It follows from Fig. 2 that, for a given frequency, Q is 
mainly sensitive to ice particles within a certain range of 
sizes. This feature is referred to as “resonance” in the fol- 
lowing discussions. The magnitude of the resonance peak 
increases with increasing frequency (f). The peak position 
in particle size, denoted by D}, hereafter, also changes with 
frequency while the ratio of D}, and wavelength à remains 
nearly constant (~ 0.4). This feature tells us that high fre- 
quencies are more sensitive to smaller particles and lower 
frequencies are more sensitive to larger particles. Especially, 
when particle characteristic size Dm is less than ~120 um, 
there exists a monotonic relationship between the polariza- 
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Circular Cylinder 7 


Polarization Difference Q (K) 
p 


0 200 400 600 800 1000 
Length of Median Mass Cylinder D, (um) 


Fig. 3. Simulated polarization differences for cirrus clouds com- 
posed of randomly oriented circular cylinders with differing aspect 
ratios (denoted close to curves) at 683 GHz. Note, the particle char- 
acteristic size Dm represents the length of the median mass cylinder 
in this case. 


tion difference Q and frequency f. Consequently, multi- 
frequency measurements of the polarization difference Q are 
potentially useful in determining ice particle characteristic 
size of cirrus clouds. However, this conclusion is limited to 
clouds with nearly spherical particles in random orientation. 


If cirrus clouds are composed of other forms of particles, 
the results will be somewhat different. In Fig. 3 simulated po- 
larization differences for circular cylinders are shown versus 
particle characteristic size. These results are for a fixed fre- 
quency of 683 GHz but with changing particle aspect ratios. 
Clearly, the aspect ratio of the ice particles plays a crucial 
role in the relations between the polarization difference Q 
and the particle characteristic size Dm. With the increase of 
particle aspect ratio, the magnitude of the resonance peak in- 
creases and the position of the peak moves toward large par- 
ticles. When D, is less than ~150 um, Q is approximately 
inversely proportional to the aspect ratio, while D,, is larger 
than ~700 um, Q becomes approximately directly propor- 
tional to the aspect ratio. Simulations for circular plates show 
similar results (see Fig. 4). However, difference is obvious 
for the two types of particles in the following two aspects. 
First, for plates the position of the resonance peak moves 
much more slowly with the increase of ice particle aspect ra- 
tio. Secondly, the magnitude of the resonance peak decreases 
slightly when the aspect ratio of plates changes from 20 to 30. 
This phenonenon does not appear for cylinders with aspect 
ratios less than 7. However, the two ranges of Dm, where a 
monotonic relationship exists for the polarization difference 
and the particle aspect ratio, are nearly identical for the two 
types of particles. 
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Polarization Difference Q (K) 
NN 


(0) 200 400 600 800 1000 
Diameter of Median Mass Plate D,, (um) 


Fig. 4. Simulated polarization differences for cirrus clouds com- 
posed of randomly oriented circular plates with differing aspect ra- 
tios (denoted close to the curves) at 683 GHz. Note, the particle 
characteristic size Dm represents the diameter of the median mass 
plate in this case. 


3.3.2 Effects of particle’s orientations and shape mixtures 


Natural ice clouds are a mixture of irregular ice particles with 
differing shapes and orientations. In this section we present 
two simple examples, in order to examine how these factors 
affect the polarization difference Q. 

In examining particle orientation effect, we follow the pro- 
cedure of Aydin and Tang (1997). Ice particles are assumed 
to be in simple forms, circular cylinder or circular plate, with 
the same aspect ratio. We assume that they have a preferred 
orientation angle relative to the horizontal plane but the par- 
ticle axis is in a uniform distribution in the azimuth. The ori- 
entation angle 6 is the angle formed by the axis of the largest 
dimension of the particle and the horizontal plane. For our 
simulations, orientation angle ô follows the truncated Gaus- 
sian function, 


1-6-5? 208? (12) 


p(s) = C 


with the normalization constant Co expressed by 
m/2 a 
Co = f eo -0/208 Sin ôdô. (13) 
0 


The particle orientation angle ô changes in the range of 
[0, 2/2]. The average orientation angle 5 is kept in zero 
in our simulations. This means, the ice particles lie in av- 
erage on the horizontal plane. When the width parameter os 
changes from zero to infinity, ice particle orientation changes 
from one preferred direction gradually to randomness. Note, 
SEP p(s) sin ddd = 1. 

Our simulations show that, when particles are oriented in 
a preferred direction, the polarization difference Q is gen- 
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Fig. 5. Effect of particle orientation: The particle orientation angle ô follows a truncated Gaussian function described by Eqs. (12) and (13). 
The average orientation angle ô is set to zero while the width parameter os changes from zero to infinity. os = 0 means fully horizontal 
orientation and og = oo means random orientation. The frequency is fixed at 683 GHz. Note, the value of Q shown here represents the 
polarization difference (in Kelvin) per 6 gm7? of IWP. Diamond marks the position of the resonance peak. 
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Fig. 6. Effect of particle shape mixture: Ice particles are a mixture of cylinders with differing aspect ratios from 1 to 7 and the number of each 
kind of cylinders is determined using a truncated Gaussian function (see details in text). The results shown correspond to width parameter oy 


changing from 0.1 to infinity. The frequency is fixed at 683 GHz. Note, Q should be taken as polarization difference (in Kelvin) per 60 gm7 
of ice water path. Diamond marks the position of the resonance peak. 


erally larger than in the case of a random orientation. Espe- 
cially, when øs is zero, which means all particles are horizon- 
tally oriented, the value of Q at the resonance peak is several 
times larger than that in a random orientation: 8 times for 
circular cylinders and over 10 times for circular plates (see 
Fig. 5a and 5b). In general, orientation effect is more sig- 
nificant for plates than for cylinders. However, with the in- 
crease of width parameter o5, Q decreases rapidly. When 
os equals 30°, Q at the resonance peak is already less than 
50% of the value for the case of øg = 0. This is true for 
cylinders as well as for plates. It is noteworthy that the peak 
position D}, of Q in particle characteristic size changes little 
for plates when the orientation width parameter os changes 


www.atmos-chem-phys.org/acp/3/39/ 


2 


from zero to infinity. For cylinders, however, D}, initially re- 
mains nearly unchanged as øs increases from zero to 20°. As 
05 approaches infinity from 20°, D}, moves gradually from 
~500 um to ~600 um. Based on this simulation, we state 
that what we found for randomly oriented particles in the res- 
onance feature of the polarization difference Q with respect 
to frequency, particle aspect ratio, and particle characteristic 
size, should remain to be similar for oriented particles. The 
only exception is that the magnitude of the resonance peak 
for oriented particles can be much larger than that of ran- 
domly oriented ice particles, depending on the width param- 
eter os. We also noticed that particle orientation has a much 
more prominent effect than particle shape on the magnitude 
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of the resonance peak (compare Figs. 5a and 5b). 

The effect of particle shape mixture is shown in Fig. 6, 
where the simulation results are given for circular cylinders. 
In this simulation, we assume ice particle numbers to be in 
a modified Gamma distribution in respect to particle size, as 
given in Eq. (11). However, ice particles at a certain size, 
denoted by the length of the circular cylinder, are composed 
of cylinders having equal lengths but with differing aspect 
ratios. The number fraction of ice particles in each aspect 
ratio is described by a discrete, truncated Gaussian function, 


g(n) = Í p=0-0)?/20;?, (14) 
Do 


n denotes aspect ratio and takes discrete values from | to 7. 
Do is a normalization constant, enabling ies 1g) = 1. 
The average aspect ratio 7 increases linearly from | to 7 as 
the cylinder length increases from zero to 2.0mm. The de- 
gree of shape mixture is controlled by the width parameter 
o,. Apparently, when oy is small, the particles at a certain 
length are composed mainly by cylinders with an aspect ratio 
corresponding to that length according to the linear relation- 
ship between 7 and particle length. When o, becomes larger 
and finally approaches infinity, particles will be composed of 
equal portions of cylinders with aspect ratios from | to 7. 
The results in Fig. 6 show that shape mixture has influ- 
ences on both the magnitude and the position of the reso- 
nance peak of the polarization difference Q. Compared with 
the results in Fig. 3, Q in Fig. 6 shows a considerably higher 
wing on the right hand side of the resonance peak, which 
extends to the large particle sizes. With an increasing width 
parameter o,,, the behavior of polarization difference Q grad- 
ually becomes close to the case of cylinders with a small but 
fixed aspect ratio (compare Fig. 3). This is understandable. 
Since the volume of a cylinder with smaller aspect ratio is 
much larger than the volume of a cylinder of the same length 
but with a larger aspect ratio, particles with small aspect ra- 
tios will gradually dominate the volume of the ice particle 
ensemble, with the increase of the width parameter o3. As a 
result, they will also gradually dominate the scattering effect. 
Although shape mixture has a considerable influence on 
the polarization difference Q, the basic feature of Q in re- 
spect to the particle characteristic size does not change. This 
means, the polarization difference at a given frequency still 
shows its resonance for particles at a certain range of sizes, 
although the magnitude and position of the resonance peak 
may change depending on how the ice particles are mixed. 


4 Conclusions 


The main purpose of this study is to promote a general un- 
derstanding of the potential of polarization measurements 
from space at mm/sub-mm wavelengths in determining cirrus 
clouds parameters but not to quantify the polarization signa- 
ture of realistic cirrus clouds. To this end, simplified radiative 
transfer procedures have been pursued, which include 
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(1) using single scattering to describe the interaction of the 
ice particles and atmospheric radiation; 


(2) ignoring the water vapor screening effect in the layers 
above the clouds; and 


(3) using circular cylinders, circular plates and nearly 
spherical particles to simulate natural ice particles. 


Another important simplification included in our simulation 
is to omit the polarization signal from the ground surface, 
which may be significant for dry atmospheres especially at 
the millimeter wave frequencies. 

Despite these simplifications introduced in the simulation, 
our results clearly show the unique features of polarization 
measurements. The most outstanding one is the resonance 
of the polarization difference for ice particles in a certain 
range of sizes. This resonance feature is also reported re- 
cently in a simulation and measurement study by Troitsky 
et al. (2001) for precipitating clouds at lower frequencies 
(<85 GHz). Furthermore, we found that the position of the 
resonance peak in respect to particle size changes from large 
to small when the observing frequency increases from low 
to high. The position of resonance peak depends also on 
the particle shape and the aspect ratio. This resonance fea- 
ture appears even when the particles are oriented or mixed 
with different shapes, although the exact magnitude and po- 
sition of the resonance peak may change accordingly. This 
unique features of polarization difference with respect to par- 
ticle microphysical properties suggest that polarization mea- 
surements from satellite at some well-chosen mm/sub-mm 
frequencies will benefit the determination of cirrus cloud pa- 
rameters. In addition, polarization measurements in mm and 
sub-mm wavelengths are also an excellent complement to 
similar measurements in infrared and visible wavelength in 
determining microphysical properties of a wide range of cir- 
rus clouds (Liou and Takano, 2002). 

The results presented in the paper have been concentrated 
on a specific looking geometry of the satellite, i.e. the satel- 
lite zenith angle was fixed at 54°. In fact, we also examined 
the zenith angle dependence of the polarization difference Q, 
although the examination was limited to randomly oriented 
ice particles. We found that, for unit IWP along the observa- 
tion direction, Q begins with a monotonic increase with the 
increase of the zenith angle. At about 60° to 70° of the zenith 
angle, Q reaches its maximum. Q goes down continuously 
from the maximum when the zenith angle approaches 90°. 
However, for a given cloud, the IWP along the observation 
path increases drastically with increasing zenith angle, there- 
fore Q generally shows a monotonic increase trend when the 
satellite zenith angle varies from 0° to 90°. 

Caution must be taken when trying to compare realistic 
measurements with the simulation results given in this study, 
especially when ice clouds are optically thick. Since multiple 
scattering has not been considered in the radiative transfer 
calculations, our simulation results are limited to optically 
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thin clouds. For example, at 683 GHz an ice cloud with an 
IWP of 6 gm~? has an optical depth of less than 0.1, when 
the cloud is composed of circular ice cylinders with an as- 
pect ratio of 5 and a median mass length of 500 um. If the 
single scattering radiative transfer can be taken to be valid 
for this case, then the polarization difference measured by a 
space-borne radiometer ranges from a few tenths of a Kelvin 
for randomly oriented particles to a few Kelvins for fully 
horizontally oriented particles (referring to Fig. 5a). As to 
the multiple scattering effect, our understanding is following. 
We should expect a linear increase of Q with the increase of 
cloud IWP in the range from zero to a certain amount, at 
which the multiple scattering effect begins to become signif- 
icant. Beyond that amount of IWP, Q should gradually ap- 
proach its saturation limit. In this process, the position of the 
resonance peak may change due to the variation of cloud op- 
tical depth with the change of ice particle size, even when the 
cloud IWP is fixed. The screening effect of the atmosphere 
in the layer on the top of the cloud is found to be negligibly 
small in the specific case which is discussed in this paper, 
since the transmittivity of this layer is high. At the highest 
frequencies (874 GHz) of our study it is larger than 0.95. 

Future work will be focused on the quantitative determi- 
nation of the polarization difference for various ice clouds 
composed of more realistic ice particles. Especially, the ef- 
fects of multiple scattering effect and cloud inhomogeneity in 
both horizontal and vertical dimensions will be taken into ac- 
count. For this purpose, a three-dimensional, fully-polarized 
radiative transfer model is under development in the Institute 
of Environmental Physics, University of Bremen (Kerridge 
et al., 2002). A cloud parameter retrieval scheme taking into 
account of the polarization measurements will also be devel- 
oped. 
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